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INTRODUCTION 


The effort during this period has concentrat’d on reduction of data 
already in hand. This includes data mainly for two projects: the abso- 

lute brightness temperature measurements at 86.1 GHz and the spectral 
line observations of Venus. 

The latter, mainly CO observations at 115 GHz, have been assembled 
into a common data base and are being searched for variation patterns. 

We plan to write up the results for submission this fall to a journal. 
The results of the former experiment have been written up and submitted 
to IEEE - Transac t Ions on A nten nas and P ropagat 1 o n . A copy of the manu- 
script is attached. 
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Abstract - Careful observations have been made at 86.1 GHz to 


derive the absolute brightness temperatures of the Sun (791A ± 

192 K), Venus (357.5 ± 13.1 K) , Jupiter (179.4 ± 4.7 K) , and 
Saturn (153.4 ± 4.8 K) with a standard error of about 3 %. This 
is a significant improvement in accuracy over previous results. 

A stable transmitter and novel superheterodyne receiver were 
constructed and used to determine the effective collecting area 
of the MWO 4.9 m antenna relative to a previously calibrated 
standard gain horn. The thermal scale was set by calibrating 
the radiometer with carefully constructed and tested hot and 
cold loads. The brightness temperatures may be used to establish 
an absolute calibration scale and to determine the antenna aper- 
ture and beam efficiencies of other radio telescopes at 3.5 mm 


wavel ength. 


I. INTRODUCTION 


The ebsol'ite intensity of a celestial radio source is difficult 
to measure accurately at millimeter wavelengths, primarily because of 
instrumental difficulties, but it conveys more information useful to 
theorists than any other observable parameter of electromng;netic 
radiation. The short wavelengths employed require extreme pre- 
cision in the fabrication of microwave hardware. The relatively 
high losses of typical waveguide components and the poor performance 
of detectors combine to significantly increase the noise temperature 
of millimeter wavelength radiometers over their lower frequency 
counterparts. In addition, the precise surfaces required for filled- 
aperture antennas at short wavelengths severely liDiit the effective 
collecting area of existing radio telescopes. The Earth's atmosphere 
will also significantly attenuate celestial radio signals due to 
absorption by molecular oxygen and water vapor {!>. The net result 
of these difficulties is a relatively imprecise absolute calibration 
scale for millimeter-wavelength radio astronomical observations. The 
purpose of this experiment is to improve the accuracy of the absolute 
flux density scale at 3.5 mm wavelength by employing new techniques 
for antenna gain calibration, thermal noise calibration, and impedance 
mismatch corrections, and by the construc*^ion of a stable transmitter 
and receiver specifically designed to reduce the dominant measurement 


errors. 


II. INSTRUMENTATION 


The system used to detect celestial radio emission consists of 
a conventional paraboloidal reflector antenna with a superheterodyne 
Dickc radiometer located at the prime focus. The thermal emission 
from blackbody loads is used to calibrate the absolute intensity of 
the detected signals. The effective collecting area of the paraboloid 
is determined by comparing the signal received from a continuous wave 
(CW) transmitter located some distance from the receiving equipment with 
the signal from a previously-calibrated standard gain horn. The same 
general procedure has been previously used to provide absolute flux 
density calibrations of radio sources at 2.3 GHz {2>, 20.5-35.5 GHz 
i3>, 35.0 GHz i4>, 97.1 GHz and lAl GHz {5>. 

A. Paraboloidal reflector 

The antenna is the 4.88 m diameter, ^/0.5 paraboloidal reflector 
<61, {7> operated by the Millimeter Wave Observatory on Mt. Locke, 

Texas at an altitude of 2015 m. In order to minimize atmospheric prop- 
agation loss, a high-altitude telescope site was chosen in an arid, 
cloud-free climatological zone. The RMS random telescope tracking 
error is 2 arc seconds, and the corrections to measured intensities 
due to tracking jitter are negligibly small (= 0.04 %) . The systematic 
pointing error is 13 arc seconds RMS, and significant intensity errors 
cun result if this pointing error is not corrected or taken into account 
properly. For these observations a five point grid centered on the 
nominal source position was mapped in order to measure the systematic 
pointing residuals and to correct the measured intensities. The re- 
flector profile has an RMS deviation from a paraboloid of only about 
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0.1 mm, insuring good aperture efficiency and nearly diffraction-limited 
performance throughout the millimeter band. The radiation patterns 
measured at 86.1 GHz (wavelength X ■ 3.A82 ram) on the pattern range 
are shown in Fig. 1. 

B. Calibration horn 

The calibration horn used as a gain standard has '>een previously 
described by Ullch {8>. It is a simple conical horn with an axial 
length of 121.92 cm and an aperture of 10.16 cm. At 86.1 GHz the 
far-field directivity was determined to be 37.202 ± 0.030 dB (all 
errors quoted in this paper are one standard deviation) using a 
radiometric calibration scheme similar to the "artificial moon" tech- 
nique f9>. This conical horn was also used to measure the average 
disk brightness temperature of the Sun and to perform sky "tipping" 
scans to monitor atmospheric extinction. Measured radiation patterns 
at 86.1 GHz are shown in Fig. 2. Fig. 3 is a plot of the spillover 
efficiency (l.e., the normalized fraction of radiated power 
which is contained within a cone of specified half-angle) calculated 
from the measured patterns. At 86.1 GHz the mnsured horn dissipation 
loss is only 2.0 ± 0.5 5!, which agrees with the calculated alue of 
2.3 Z {8>. 

C. Transmitte r 

A CW transmitter was constructed to provide a strong signal for 
the relative directivity measurements of the paraboloid and the cali- 
bration horn. A phase-locked klystron at 86.100 GHz provided power to a 
lens-corrected conical horn antenna with about 33 dB gain. As shown 
in the block diagram of Fig. 4, a feedback control circuit was used 
to maintain constant transmitted power. Laboratory mo.'isurcments 


indicate a typical stability of about 0.04 dB/hour. The entire 
transmitter is contained in a thermally insulated box to reduce errors 
caused by temperature fluctuations. The transmitting horn was carefully 
aligned to the direction of the telescope to insure uniform illumination 
of the paraboloid. 

D. Receiver 

A multipurpose receiver was also constructed to enable three types 
of observations to be conveniently made. Fig. 5 is a block diagram 
showing the significant subsections of the receiver. The first mode of 
operation is that of a single-conversion, phase- locked, narrowband super- 
heterodyne receiver with an intermediate frequency (IF) of 30 MHz. An 
accurate waveguide below cutoff (WBCO) attenuator (Airborne Instruments 
Laboratory Type 3232) is used with the IF substitution technique fl0> to 
determine the ratio of the transmitter power received by the paraboloid 
to that received by the standard gain horn. The accuracy of the WBCO 
attenuator calibration is ±{0.03 dB + (0.005 dB per 10 dB change in reading)!. 
The waveguide switch at the receiver input allows rapid comparison of 
the signals from the paraboloid feed horn and from the calibration horn 
which is physically coaxial with the telescope axis. The WBCO attenuator 
is preceded by a high-power, heat-sunk transistor amplifier (+29 dBm for 
-1 dB compression) to insure adequate djTiamic range to acconmiodate the 
33 dB difference in received signal lev^l without significant IF satur- 
ation. A 30 ^^^z bandpass filter reduces the bandwidth of the detected 
spectrum and Improves the signal-to-noise ratio. In operation, a change 
in radio frequency (RF) signal level is compensated by UTiCO attenuation 
of the IF signal to maintain a constant detector output. All waveguide 
flanges and other sites of possible RF leaV.age are wrapped with copper 
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foil tape to insure sufficient isolation between the paraboloid and the 
calibration horn. Special observations were made in an effort to detect 
residual leakage, but none was found at a level -23 dB below the cali- 
bration horn signal. 

The second mode of receiver operation is that of a broadband 
Dicke superheterodyne radiometer. In this mode the IF bandpass is 
from 10 to 500 MHz. The receiver noise temperature was measured to 
be 6100 K (double sideband) at the input, and the was 0.6 ± 0.1 K 

for one second of integration time, which agrees with the theoretical 
value of 0.55 K. A switchable circulator alternately connects the 
receiver to one of two identical feed horns at a 20 Hz rate. The feeds 
are optimized pyramidal horns which have aperture dimensions of 0.90 X 
in the ^-plane and 1.25 X in the H^-plane. The feed spacing resulted 
in an orthogonally polarized reference beam separation of 2.6 half- 
power beamwdiths off-axis in hour angle. Observations of blackbody 
calibration loads and of celestial radio sources were made using the 
broadband radiometer output. 

The third mode of receiver operation is that of a ref lectometer . 

An integral waveguide slotted line was used to measure the reflection 
coefficients of the various antennas, of the thermal calibration loads, 
and of the receiver input. 

E. Thermal calibration loads 

Absolute noise power spectral density calibration was obtained 
by observing the difference In signal level from two loads. Both 
loads were fabricated using Eccosorb CV-3 porous carbon- impregnated 
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foam absorber (Lmerson and Cuming, Inc., Canton, MA) enclosed in an 

aluminum cavity -Clll. One load was nominally at ambient temperature, 

and a laboratory thermometer inserted in the absorber was used to 

determine its physical temperature. The second load was Insulated 

on all sides with 5.1 cm of expanded polystyrene and filled with 

-3 

liquid nitrogen. The foam walls (density “ 0.032 g cm ) have an 
estimated relative dielectric constant of 1.029 and a loss tangent of 
about 3.6 X 10 ^ near 3 mm wavelength based on measurements of poly- 
styrene at 90 GHz {12> and density scaling laws <13>. The calculated 
power reflection of the foam wall is 0.02 X and the absorption is 
0.33%. The measured transmissivity was 0.9969 ± 0.0007 at 86.1 GHz, 
which agrees vrithli. the error with the calculated value of 0.9965. 


111. \Sm MEASUK»IENTS 


The waveguide slotted line built into the receiver was used to 
measure the voltage standing wave i*atlo (VSTO) of system components, and 
the results are given in Table 1. The magnitudes of the voltage re- 
flection coefficients (Ir]) at 86.1 GHz are generally quite low. The 
receiver input is the worst match of the components tested even though 
an isolator (see Fig. 5) has been used to improve the impedance matching. 
Measurements of the mixer without the Isolator indicated an input VSWR of 
2.7, 80 it appears that the isolator dees significantly improve the 
matching. Tlie reflection coefficient of the /i.9 m paraboloid was de- 
termined by axially focusing the receiver box. Theoretical calculations 
{1A> predict a value of 0.0012, which agrees with the measurement. The 


Source 


Calibration horn 
Reflector feed horn 
A. 9 M Paraboloidal reflector 
Ambient temperature load 
Eccosorb CV-3 absorber 
Liquid nitrogen cooled load 
Receiver waveguide switch input 
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ILE I 

VSWR DATA 


VSWR 


0.029 

4 

0.001 

1.060 

± 

0.002 

0.024 

4 

0.001 

1.049 

4 

0.002 

0.0015 

4 

0.0C05 

1.003 

± 

0.001 

0.0020 

4 

0.0005 

1.004 

4 

0.001 

0.0020 

4 

0.0005 

1.004 

1 

0.001 

0.044 

4 

0.002 

1.092 

4 

0.004 

0.076 

± 

0.003 

1.165 

4 

0.006 
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thermal calibration loads and the absorber were measured using a 
traveling mechanical stage In front of the calibration horn. These 
data were taken at normal incidence and corrections due to waveguide 
loss were determined by observations of a copper shorting plate. 
Measurements of the liquid nitrogen cooled load with cryogen 
(|r| ■ 0.044 ± 0.002) and without cryogen (|rl ■ 0.011 ± 0.002) 
indicate that the dominant source of reflection is the foam-cryogen 
interface. 

All pattern range, thermal calibration, and radio source measure- 
ments have been corrected for mismatch loss -CIS}. However, since only 
the amplitudes of the reflection coefficients (and not the phases) arc 
known, an "exact" correction cannot be made in general. We have 
chosen to set the conjugate mismatch loss factor equal to the mean of 
its maximum and its minimum possible values and, for the purpose of 
conservative error analysis, to estimate an error of one standard 
deviation by the maximum possible error. Thus the mean mismatch loss 
correction factor between components 1 and 2 is (assuming small re- 
flection coefficients) 


<M> 


Power dissipated in load ^ 

Maximum available power ^ 1 ^ t \ ^ > 


and the estimated lo error is 



( 1 ) 


( 2 ) 


r 
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The effect of mismatch on the broadband radiometric data (instan- 
tai. ms bandwidth 1 CHz) is less seric s than for the coherent 
measurements since the detector will in n 'st cases average the signal 
over several cycles of the sinusoidal standing wave pattern, and the 
mismatch '"actor given by (1) will then be nearly correct. In this 
case the ' rue mismatch error will be overestimated by (2). 


IV. PATTERN RANGE MEASUREMENTS 

An unobstructed pattern range was used to provide a quasi-plane 

wave input to the 4.9 m paraboloid and to the coaxial calibration 

horn for the narrowband directivity comparison measurements. The 

transmitter site was a hilltop in the Davis Mountains State Park at 

an altitude of 1680 m and a range of 12.9 km from the MWO. As seen 

from the telescope, the transmitter is at -1.5° elevation angle, and 

the closest possible ground reflection is about 0.63° below that angle 

(i.e., about 13 beamwidths off-axis). A search for signals near this 

angle indicated that ground reflections w’ere negligibly small (less 

than -30 dB with respect to the direct path). The finite range of the 

transmitter, however, means that the wave impinging on the telescope 

is not planar, but rather spherical, and the resulting quadratic phase 

error produces some reduction in peak gain of the paraboloid. Fortunately, 

a degree of freedom exists in the telescope which can be utilized to 

simulate far-zone radiation patterns in the Fresnel zone. At a range 

2 

R = 12. km the transmitter is 1.89 D /A from the paraboloid (of diameter 
D). By defocusing the feed horn axially av/ay from the reflector (of 
focal length f) by a distance c, the far-zone radiation pattern can be 
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2 

virtually duplicated for ranges greater than 0.5D /X. The axial de- 
focusing distance c is given by 


c 



1 

4(f/D) 



(3) 


where K ■ 0.92 {161 . The peak gain of the antenna is the same whether 
it is focused in the Fresnel zone or at infinity {171 . First the feeds 
were focojed to receive maxiiEum transmitter signal and the directivity 
comparison data were taken. Then the feeds were axially displaced by 
0.15 X (calculated from (3)) toward the reflector to focus at infinity, 
and the planetary radio observations were made at this setting. In this 
case the expected gain reduction if the feeds were not refocused is only 
about l.A Z. 

Figure 6 is a plot of the variation in signal power received by 
the paraboloid during a typical 20 minute period. Most of the noise is 
due to atmospheric scintillation, and data were taken between midnight 
and sunrise to niniml;;e errors due to fading. Table II presents individ- 
ual measurements of the obsers’ed signal power ratio. Some of these 
measurCkTkentw were made at lower transmitter power levels to search for 
possible receiver nonlinearity due to saturation, but in no case was 
any nonlinearity detected. A calibrated rotary-vane RF attenuator in 
the transmitter was also used to check the WBCO IF attenuator calibration, 
and no nonlinearity was detected within the accuracy af the RF attenuator 
calibration (- 0.1 dB) Each measurement was made after careful pointing 
adjustments to maximize the received signals. 
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TABLE II 

ANTENNA SIGNAL RATIO DATA 


Trial 

Date 

UTC 

(Hours) 

SlRiinl from reflector 
Signal from horn 
(dB) 

1 

12/14/77 

0500 

32.914 

2 

12/14/77 

0900 

32.885 

3 

12/15/77 

0400 

32.924 

• 

12/15/77 

0900 

32.928 

5 

12/15/77 

1400 

32.930 

<l-5> 



32.916 i 0.008 




lA 


The reflector directivity (referred to the aperture plane of the 
on-axis feed) is simply obtained by correcting the observed signal 
power ratio for the different losses of the paths from the feed horn 
and the calibration horn to the ferrite switch and then multiplying 
by the directivity of the calibration horn. As Indicated in Table III, 
the directivity of the A. 9 m reflector at 86.1 GHz is 70.160 ± 0.072 dB. 
The estimated uncertainty is the quadrature sum (root sum of squares) 

of the estimated errors (which are assumed to be independent) in the 

$ 

calibration horn directivity, signal power ratio repeatability, mis- 
match correction, receiver gain correction repeatability, WBCO atten- 
uator calibration, polarization angle misalignment, refocusing, anJ 
difference in cross polarization. This last uncertainty exists be- 
cause the calibration horn directivity was determined using an un- 
polarized source, whereas for the directivity ratio experiment a 
linearly polarized source was used. Table IV summarizes the measured 
parameters of the calibration horn and of the paraboloid. The reflector 
has a high aperture efficiency (53.58 ± 0.89 %) , indicating that 
surface errors are much smaller than a wavelength and that the radia- 
tion pattern is nearly diffraction-limited. The reflector half power 
beamwidths are also essentially equal since the aperture dimensions of 
the pyramidal feed horn have been chosen to produce a imain lobe with 
a circular cross section. 


V. THERMAL CALIBRATION 


The noise power P (W) delivered to the receiver in a single polarization 


through the aperture of either the on-axis feed horn or tlie calibration 
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TADLE III 

REFLECTOR DIRECTIVITY 




Quantity 

Value 

(dB) 


f 


i 


Calibration horn directivity 

Ratio of reflector signal to horn signal 

Uncertainty due to mismatch 

Ratio of receiver gain at calibration horn 
aperture to gain at reflector feed aperture 

Uncertainty in WBCO attenuation 

Uncertainty due to polarization angle mis. ' ignment 
Uncertainty in refocusing 

Uncertainty due to difference in cross polarization 
Reflector directivity at horizon 


37.202 -t 0.''10 
32.916 i 0.008 
± 0.026 

0.04 1 0.018 


± 0.046 
i 0.005 
± 0.025 


± 0.022 


) 


70.160 ± 0.072 


I 

I 


I 


f 

L 

^ ' 


TABLE IV 


ANTENNA PARAMETERS AT 86.1 GHZ 


Quantity 


Calibration Horn 


Diameter (m) 

2 

Geometrical area (m ) 
Directivity (dB) 

2 

Effective area (m ) 

Aperture efficiency (2) 

HPBW in E-plane (arc minute) 
HPBW in H-plane (arc minute) 


0.1016 
0.008107 
37.202 ± 0.030 
0.005065 ± 0.000030 
62.48 ± 0.44 
127 ± 2 


Paraboloidal Reflector 


4.677 

18.679 

70.160 ± 0.072 

1 

10.010 ± 0.166 
53.58 i 0.89 
2.91 ± 0.03 


152 ± 2 


2.86 ± 0.03 
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horn from a perfectly outtched (blackbody) load at an abaolute physical 
temperature T enclosed in a perfectly conducting cavity is given by 

P - k B J(T) (A) 

-23 -1 

where k is Boltzmann's constant (1.38034 x 10 J K ) , B is the 
predetection bandwidth (Hz), and the effective radiation temperature 
J(T) is given by 

J(T) ■ 

cxp(hv/kT) -1 
-34 

In ') h is Planck's constant (6.6256 x 10 J s) and v is the frequency 

(Hz). In practice tJie bandwidth B is usually not precisely known, and 
microwave radiometers are generally calibrated in terms of the effective 
radiation temperature J(T), which is proportional to iv'ise power spectral 
density. When hv<<kT, J(T) “ T, and the familiar concept of antenna 
temperature may be used in the Kayleiph-Jeans approxim.it ion to the Planck 
distribution. 

The broadband radiometric observations of the Sun and planets were 
calibrated by first comparing the signal from a neon noise tube to that 
from the hot and cold loads. Tlie calibrated noise tube signal was then 
used as a secondary transfer standard to allow rapid and convenient cal- 
ibration of celestial radio sources. In practice the noise tube was 
constantly on to avoid thermal output drift, and a ferrite modulator 
was turned on and off under computer control to generate the calibration 
signal. Table V presents a summary of the noise tube calibration measure- 
ments, which were interleaved in time with the radio source observations 
to reduce systcm.it ic errors. Alternate observations of the hot and cold 
loads over the apertures of the on-axis feed horn and of the calibration 
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horn were also used to dctcrnlne the relative receiver gain from these 
two input ports to the square law detector. As indicated in Table V, 
simultaneous recordings of ambient (hot load) temperature T.„_ and 
of barometric pressure were used to determine the radiation temperature 
difference of the two loads. The boiling point of liquid nitrogen 
*^LN related to the barometric pressure (mm Hg) by {18> 

and the microwave absorber is physically at this temperature since it 
is completely immersed in the cryogen. The foam wall, however, 
attenuates the emission from the absorbei and also emits radiation 
according to its normalized power loss factor a. For small lossr 
the effective radiation temperature J (T^) of the cold load is 

J(Tc> - + u/2) (7) 

where a linear temperature gradient in the foam wall has been assumed 
and reflections are neglected <18>. In the present case a - 0.0033 
and the foam '.’all increases the effective cold loud temperature by about 
0.3 K. This was experimentally verified by comparison measurements of 
a second cold load without the foam wall. A A3° flat reflector was 
placed above nn absorber-lined cavity filled with cryogen, and the 
signal was found to be within 1 0.4 K of the cold load with foam. 
Measurements made with a second foam wall in place indicated an in- 
crease of 0.9 1 0.4 K, which agrees with the calculated value of 0.7 K. 
Thun the low loss of the foam wall and its small effect on the cold 
load temperature has been confirmed. 


12/15/77 

12/17/77 

12/18/77 

12/19/77 

12/19/77 


TABLE V 

THERMAL CALIBRATION DATA 



Barometric 


UTC 

(Hours) 

Pressure 
(^Ci Hg) 

Power Output 
(Volt) 

T 

AMB 

(K) 

^LN 

(K) 

0000 

609 

0.923 

281.3 

75.7 

0200 

603 

0.917 

28^4.4 

75.6 

0300 

605 

0.920 

280.1 

75.7 

0300 

609 

0.938 

280.5 

75.7 

0000 

606 

0.928 

284.0 

75.7 

1600 

607 

0.925 

284 . 4 

75.7 


^AMB~ ^LN 
(K) 

&CAL 

WC (Feed) 

AH/C (Feed) 
AH/C (Horn) 

^CAL 

(K) 

205.6 

0.1389 

1.0055 

28.30 

208.8 

0.1428 

0.9945 

29.53 

204.4 

0.1459 

0.9933 

29.54 

204.8 

0.1426 

0.9760 

28.94 

208.3 

0.1393 

0.9840 

28.76 

208.7 

0.1405 

0.9889 

29.05 

- 

- 

0.9904 ± 0.0041 

29.02 ± 4 ^ 
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Mismatch will also increase the effective temperature of the 
cold load since some noise power radiated by the receiver out its 
input port will be reflected back into the receiver. The presence 
of an Isolator in the receiver input waveguide assures that the 
noise power radiated out the receiver will be at nearly ambient 
temperature, and this was confirmed by observations of a shorting 
plate at the feed aperture. The effect of mismatch between the 
cold load and the receiver input (evaluated using (1) and adding the 
reflected signal) is to increase the effective radiation temperature 
of the cold load by about 1.6 K, while having essentially no effect 
on the ambient temperature load signal. Thus the cold load radiation 
temperature is approximately equal to J(T..,) + 1.9 K ■ J(77.8 K) fo’' 
typical conditions on Mt. Locke. The average value of the effective 
noise tube signal from Table V is 29.02 ± 0.19 K referred to the 
aperture of the main feed, where the error is estimated from the re- 
peatability of the measurements. In addition to this random error, 
uncertainty in the load tempcrnturcs (conservatively estimated at 
1 0.5 K each), in the foam wall loss (± 0.07 %), and in the mismatch 
correction (i 0.76 X) combine in quadrature to produce a total uncertain- 
ty in the thermal calibration scale of 1.06 X. 

VI , ATMOSPHERIC EXTINCTION 

The absorption of radio waves at 86.1 GHz by the terrestrial atmos- 
phere is both significant and variable. Therefore accurate measurements 
of celestial sources must be accompanied by a nearly simultaneous deter- 
mination of atmospheric opacity. In principle this can be done by 
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observing the decrease in received signal Intensity as the effective 
air mass of the source increases. However, in this case atoiospheric 
absorption cannot be separated fron changing telescope gain (due to 
gravitational distortion). Furthermore, significant changes in atmos- 
pheric absorption can occur during the several hours usually required 
for such observations. Measurements of atmospheric emission, however, 
can be used to directly determine atmospheric absorption due to oxygen 
and water vapor. The two-layer model of Kutner {19] was used with 
ground temperature data to calculate the effective radiation temperatures 
of oxygen and water vapor. Surface absolute humidity measurements 
were also made to estimate the (variable) precipitablc water vapor 
in a column and thus the relative absorption of water vapor {20> com- 
pared to the essentially constant oxygen component. From these values 
a weighted mean atmospheric radiation temperature was calculated and 
recorded for future use. The actual sky brightness temperature data 
were taken with the calibration horn rather than with the paraboloidal 
antenna, since the horn sidclohe levels were more accurately known 
than those of the paraboloid. Of c.uirse, the signal received by the 
horn is actually a convolution over steradians of the horn power 
pattern and the brightness distribution of its surroundings. A computer 
model was constructed which included the effects of atmospheric emission, 
ground emission, cosmic microwave background emission, and reflection 
of horn backlobes by the paraboloid and by the receiver box. The 

apparent sky temperature J(T , ) is given by 

SKY 
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■’"sky* ■ •’<V <*-*■'*> ♦ <®> 

where is the calculated mean atmospheric radiation temperature {21> 
and is the microwave background brightness temperature (2.8 K) 

{22>. The zenith optical depth is T, and the air mass A is given 
(to a good approximation for A<10) by 

A - 1 ( 9 ) 

sin (L) 

where E is the elevation angle. In this model the ground emission 
temperature is assuned to be the same as the surface ambient temperature. 
The measured iadiation patterns of the calibration horn given in Fig. 2 
were used for the computerized convolution, and the values of t for 
each observing session were found by the method of least squares. 

The standard deviation of the differences between the computer model 
and the observed effective temperatures w.is typically only 0.3 K, 
indicating that the model quite accurately fits the shape of the sky 
tipping curve. Figure 7 shows the results of one such fit, and the 
values of zenith optical depth arc given in Table VI. The values of 
T are quite low (0. 052-0. 07.'») and are consistent with only 2-4 mm 
precipltuble water vapor. The error in the atmospheric absorption 
correction to observed celestial source intensities was conservatively 
assumed to be 20X of the actual correction. 

VII. DEPENDENCE OF REFLECTOR GAIN ON ELEVATION ANGLE 

A scries of observations were made on 24-27 November 1978 to 
determine if the peak gain of the paraboloid varies with elevation angle. 


TABLE VI 


ASTRONOMICAL DATA 





Date UTC Source 
(Hours) 


Nur..ber 
of Scans 


T±lo 

(Nepers) 


Antenna 


T^±lo 

(K) 


S±lo 

(Jy) 


(K) 


12/17/77 

2000 

Sun 

5 

0.075 ± 

0.015 

Horn 

229.9 

± 8.8 

(125.3 ± 4.8)xl0^ 

7985 ± 305 

12/13/77 

20C0 

Sun 

11 

0.053 ± 

0.011 

Horn 

227.0 

± 5.6 

(123.8 ± 3.0)xl0^ 

7885 ± 193 

12/18/77 

0830 

Jupiter 

10 

0.076 ± 

0.015 

Reflector 

5.572 

♦ 0.102 

1482 ± 28 

175.3 ± 3.3 

12/19/77 

0730 

Jupiter 

10 

0.052 ± 

0.010 

Reflector 

5.542 

± 0.067 

1529 ± 19 

180.9 t 2.2 

12/19/77 

1300 

Saturn 

4 

0.052 ± 

0.010 

Reflector 

0.747 

± 0.015 

206 ± 4 

153.4 ♦ 3.0 

11/25,26, 

27/78 

1200 

Jupiter 

108 

0.05 to 

0.12 

Reflector 

4.205 

± 0.061 

1160 ± 17 

180.0 ♦ 2.6 

11/24.25, 

26.27/78 

1700 

Venus 

56 

0.05 to 

0.12 

Reflector 

15.490 

± 0.420 

4273 ± 116 

357.5 ± 9.7 


W 
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In order to accurately correct for rapid variations in atnoapheric 
absorption, a rotating chopper wheal was used to produce a calibration 
signal according to the method first described by Davis and Vsud^n Bout 
{23> and later in detail by Ullch and Haas -C24}. The resulting in- 
tensity scale is "corrected antenna temperature" (T* ), and is simply 
the antenna temperature that would be measured outside the terrestrial 
atmosphere with a loasless antenna. For a source which completely 
fills the forward beam, the corrected antenna temperature equals the 
source brigtitness temperature. Since the Sun is much larger in angular 
extent than either the n^in beam or the ncor sidelobeo of the paraboloid 
power pattern, gravitational (or thcvoial) distortion of the reflector 
cannot affect the received signal by a significant amount. Thus 
observations of the Sun were used to check the accuracy of the atmos- 
pheric absorption corrections with the chopper-wheel calibration tech- 
nique. An shovm in Fig. 8, there is no indication of any change in 
measured solar intensity with elevation angle. The solar brightness 
temperature at 86.1 GHz measured with the 4.9 m poraboloid and the 
chopper-wheel calibration scheme is 7750 i 150 K using a calculated 
beam coupling efficiency of 0.77; the error is due solely to receiver 
noise and atmospheric absorption corrections. Taking into account the 
uncertainty in telescope losses, the solar brightness temperature is 
7750 i 350 K. Concurrent observations of the planets Venus snd Jupiter 
do show a slight increase In intensity from the horl;',on to the zenith, 
us shown in Fig. 9. The ratio of Venus to Jiipilcr was precisely measured 
when they were at equal elevation angles, and this ratio was used to 
place the Venus and Jupiter intensities on n common scale. Formally, 
a least squares fit of the relative antenna gain function ol the 


25 


(assumed) form 

G/G - 1 + AG sin (E) (10) 

o 

to the Sun data in Fig. 8 resulted in AG ■ 0.000 ± 0.011. For the 
planetary data in Fig. 9, AG ■ 0.020 i 0.015, which indicates a 
slightly significant (2.0 i 1.5 %) Increase in gain at the zenith 
compared to the horizon, where the pattern range measurements were 
used to determine the absolute telescope gain. 

Two other methods were used to check for gain variations v/ith 
elevation angle. The first method was to compare the antenna half 
power beamwidths as measured on the pattern range with those inferred 
from Jupiter data at about 75° elevation angle. Both and H-plane 
beamwldths were narrower (after correcting for the slight broadening 
due to the finite disk size of Jupiter) at 75° elcvatio., angle than 
at the horizon, Inferring that AG * 0.018 ± 0.018. To derive this 
result, we have assumed that gravitational distortions are comparable 
in correlation length to the diametCT- of the paraboloidal reflector 
and thus that the power received by the main lobe remains constant. 

The second method was to determine the variation in the signal from 
Jupiter observed on 18 December 1977 after correction for atmospheric 
absorption utilizing the zenith optical depth obtained from the sky 
tipping procedure with the conical horn. This method gave oG = 0.027 
± 0.023. Thus all three methods of analysis are compatible with a total 
relative gain variation AG 0.021 ± 0.010, and this weighted average 
was actually used to refer planetary intensities measured at positive 
elevation angles to the hor' on, where the absolute gain calibration 
was previously accomplished. 


VI II. ASTRONOMICAL OBSERVATIONS 
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A . lackbody disk source of temperature which is small 

in angular size compared to the receiving antenna bcamwidth will 

-2 -1 

produce a measured incremental flux density S (W m Hz ) at the 

M 

surface of the Earth given by {24} 

2kv^n C e"’* 
c 

where c Is the velocity of electromagnetic wave propagation in a 

8 —1 

vacuum (2.99793 x 10 ms ), 0^ (sr) , is the solid angle of the disk 

source, C is a normalized size correction factor to account for partial 
o 

resolution of the disk by the antenna beam, and the ot^''r parameters 
have been previously defined. The solid angle correction factor C 
is given to a good approximation (when X < 0.4) by {24} 


C 


S 



( 12 ) 


where X = (4 ln2)^ R/6, 


(13) 


and R (rad) is the geometrical mean of the polar and equatorial semi- 
diameters of the (slightly elliptical) disk source and 6^^ (rad) is the 
antenna hal^ power beamwidth. The solid angle of the source is given by 




2n { l-cos(R) } = IT R^. 


04) 


The source flux density S that would be measured with no atmospheric 
loss is given by i24> 
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S $ S„ e 


,A .iijA 

Vg 


( 15 ) 


where T is the source antenna temperature (K), n. the antenna 
A A 

aperture efficiency, and A is the antenna geometrical collecting 

2 

area (m ). For cosmic radio sources the usual unit of flux density 

•■26 “2 "1 

measurements is the Jansky (1 £ 10 W m Hz 

Observations of the Sun with the conical horn were made on two 
days in December 1977, and observations of Venus, Jupiter and Saturn 
were made with the paraboloidal reflector on 6 days in December 1977 
and in November 1978. The peak source intensities were measured 
relative to the noise tube signal which was absolutely calibrated 
against the hot and cold loads. All measurements vjere corrected for 
pointing error by fitting a two dimensional gausslan function to the 
measured five point grid map to determine the peak Intensity that 
would have heem measured with no pointing error. Saturn was observed 
alternately with the on-axls and the off-axis paraboloidal feeds to 
improve the signal-co-noise ratio. Separate measurements of Jupiter 
in each beam indicated a relative gain of the reference beam of 0.922 
± 0.010 times the gain of the main beam, and this value was used to 
normalize the Saturn observations. Calculations predict a gain re- 
duction of about -0.4 dB for a feed 2.6 bcarawidths off-axis {251, in 
agreement with the measurement. Table VI lists the antenna temperature, 
the flux density, and the brightness temperature for each data set. 

The quoted errors are the quadrature sums of the Qo) random errors 
only and do not include the correlated syst''matic errors in antenna 
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t;ain or in thermal calibration. The angular r-ourco semidiumeturs 
were calculated for the time of each observation by dividing the 
knovrn unit semidiaineters <261 by the source distance in astronomical 
units. The Saturn brightness temperature has been calculated ussuroing 
that the rings do pr.^. emit or scatter radiation frota the planet. 

Thus only the disk solid angle has been used in (14), although it has 
been corrected for polar tilt as seen from the Earth. At the time 
of these observations, the Saturnicentrlc latitude of the Earth (B) 
was -9.9^, and the calculated disk solid angle has been corrected 
slightly to allow for this viewing angle. 


IX. A Nnj FLUX DEi’SlTV SCALE 


The average sol..r and planetary brightness temperatures and 
theii total (random plus systematic) errors listed in Table Vll may 
be used with (11) to establish an absolute flux density calibration 
scale near 8b. 1 CHz. This should be pat ticul.nrly useful for cali- 
brating other millimeter wavelength antennar. for radio astronomy 
purposes. An accurate and Internally consistent intensity scale is 
needed for proper theoretical interpretation of astronomical observa- 
tions, and the purpose of this experiment was to imilic such a scale 
routinely available to all millimeter wavelength observatories. The 
accuracy of the now 3..S inm ealihratJini sc.ale may he Judged by comparing 
these m.-asuxementb with previously reported observations at nearby 
frequencies and with theoretlc.al models. The 86.1 (diz solar brightness 
temperature measured with the conical horn J.s 7914 ± 192 R, which 
agrec'S with th? value of 7750 i 350 K deliTtiined by the chopi'er-wheel 
calibration technique using the 4.9 m reflector. Llnsky {27} predicts 


TABLE VII 

86.1 GHZ DISK BRIGHTNESS TE.» 


Disk Brightness 

Source Temperature 

(K) 

Sun 7914 

Venus 357.5 

Jupiter 179.4 

Saturn (B ■ -9.9*^) 


153.4 



Total lo Error 
(K) 


± 192 

± 13.1 

1 4.7 

± 4.8 
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a brlghtnesii temperature of about 7800 ± 300 K at 86.1 CHz from his 
linear fit to a recalibrated solar spectrum. Tlius the present measure- 
ments are in good agreement with previous observations. It should be 
noted that the solar brightness temperature determined by the horn 
measurements is an average over the unresolved visible disk, and any 
limb brightening will tend to increase the disk average temperature 
over the central brightness temperature which is sampled by high 
resolution antennas. 

The 86.1 GHz brightness temperature of Saturn at B ■ -9.9° re- 
ferred to the atmospheric disk only is 153.4 K with a random error of 
3.0 K and a total error of 4.8 K. This result is in good agreement 
with both previous measurements {4>, {28> and with the theoretical 
spectrum of atm.>spheric emission -C2S}. However, the ch.'’nglng tilt 
angle of the rings as seen from the E;irth introduces a slowlyr-varying 
component to Saturn's total 3 mm emisr.ion which may be as large as 
15 K f28>. Thus Saturn is not a good calibration source for long time- 
base expel iments at n few millimeters wavelength. 

Jupiter's thermal brightness temperature is not known to be 
time-variable, and its constancy and large flux density make it the 
primary calibration source at millimeter wavelengths. Our observa- 
tions indicate a brightness temperature for Jupiter of 179.4 K at 86.1 
GHz, with fl random error of 1.5 K and a total error of 4.7 K. A 
comparison ”ith both the previously observed miciowave spectrum fAl, 
f29J and with a convective model atmosphere •t29> again shows agreement 


within the quoted errors. 
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Ac 86.1 GHz Che measured brlghcncss CctapernCure of Venus is 
357.3 K with a random error of 9.7 K and a Cocal error of 13.1 K. 

Ac Che Cime of chese observaCions the planeCocenCric phase angle 
beCween Che Sun and Che EarCh was ■ 215°. Previous observaclons 
of Venus aC 97.1 GHz {4} and ac 88 GHz <30> arc in agreement wich 
this more precise result. At 88 GHz the reported mean temperature 
of Venus must be scaled upward by Che ratio of Che observed solar 
brightness temperature (7914 ± 192 K) to the assumed solar bright- 
ness temperature (6600 K) , resulting in a Venus temperature of about 
355 1 32 K. Janssen ‘C31> has calculated the millimeter-wavelength 
spectrum of Venus based on a model atmosphere which is consistent 
with all spacecraft and ground-based microwave observaCions. At 
86.1 GHz the predicted disk brightness temperature is 363 ± 10 K, 
which also agrees with the observed temperature of 357.5 i 13.1 K. 

The measurements listed in Table VII define an absolute flux 
density scale which may be used to calibrate 3 mm radio astronomical 
observations to an accuracy (lo) of about 3 X. The measured solar 
and planetary brightness temperatures are sliown to agree with previous 
measurements at nearby wavelengths and with theoretical spectra. 
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FIGURE CAHIONS 


Fig. 1. Mcanured radiation pattoms of 4.9 m paraboloidal antenna 
at 86.1 GHz. The asynmotry of the flrat aldelobcs indicates the 
presence of coma due to a slight lateral offset of the feed horn 
from the axis of tha paraboloid, (a) In hour angle (li-plane). 

(b) In declination (E-plane). 

Fig. 2. Measured radiation patterns of conical horn at 86.1 GHz. 

(a) In hour angle (H-planc). (b) In declination (E-planc). 

Fig. 3. Spillover efficiency of conicel horn at 86.1 GHz calculated 
from measured patterns in Fig. 2. 

F-'g. 4. Block diagram of Lrensmitter with feedback circuits to 
stabilize frequency and output power. 

Fig. 5. Block diagram of receiver with narrowband and wideband IF 
systems and integral reflectometer . 

Fig. 6. Variation in transmit Ler power received by the 4.9 m antenna 
over a 12.9 k>u pattern range (1235 UTC on 15 December 1977). 
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rig. 7. Measured effrctlve a.nenna tcaperature of ihe calibraclon 
horn as a function of air nasa (1530 UTC or. 19 Deceuber 1977). 

The solid line la the computer convolution model fit to the obaer- 
Vbtiona with a zenith optical depth of 0.052 Np. With no atmosphere 
present (air mass •• 0), the effective antenna temperature of the 
2.8 K cosmic background radiation is about 1.2 K at 86.1 GHz. 

Fig. 8. Measured corrected antenna tempt raturea of the .Sun at 86.1 GHz 
with the 4.9 b paraboloid aa a function of elevation angle. The 
solid line is the result of the least squares fit of the variable 
gain model described in the text with AG •• 0.000 1 0.011. 

Fig. 9. Measured planetary corrected antenna temperatures at 86.1 GHz 
with the 4.9m paraboloid as a function of elevation angle. The 
solid line lo the result of the least squares fit of the variable 
gain model described in the text with AG 0.020 ± 0.015. 
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